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Erratum on the paper: "Influence of a Feshbach resonance on the photoassociation of LiCs", New Journal of Physics, vol. 

In the paper by Deiglmayr et al [I J, we found a sign error in the transformation 
of the scattering wave function between the atomic basis (eq. 6) and the molecular 
basis (eq. 7) that we used for the case rrip = 5 for the following atomic state: 



In consequence, the results of the left column of Fig. 5 of the quoted paper are 
correct, while those of the right column are not. The corrected hgure is now displayed 
in Figure Q] below. In particular there is not anymore influence of the triplet part of 
the wave function on the singlet one (see upper panel of the right column). This affects 
also Figures 6 and 7 of reference [T] . The black trace of the upper part of Figure 6 is 
wrong, as well as the results of Figure 7. 

The main consequence concerns the interpretation of the experimental results 
reported in Figure 8 of Deiglmayr et al [T]. The magnitude and modulation of the 
photoassociation (PA) rate observed in the experiment can be interpreted only if the 
triplet component of the scattering wave function is involved, as it was discussed in the 
paper. However, the model employed in reference [T] cannot reproduce these results. 

The role of the triplet component of the scattering wave function can be invoked 
in the PA towards the B 1 !! levels by considering the small spin-orbit coupling of the 
B 1 !! state with the neighboring triplet states (see Figure 4 of Deiglmayr et al pQ), 
while the modulations of the PA rate may be influenced by a slight variation of the 
scattering length of the lowest a 3 S + state. These issues will be considered in a future 
paper. 

We would like to thank Prof. E. Tiemann for pointing out to us this error sign. 
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Figure 1. Components of the energy-normalized radial wave function for a LiCs 
pair colliding at 500 fiK, in the atomic basis (left column) and in the molecular 
basis (right column). The displayed case corresponds to m/ 1 +m/ 2 =5. Under the 
given experimental conditions, the upper-most atomic state fi = 2, m,f 1 = 2, 
/2 = 3, mj 2 = 3 is the only open entrance channel, the three other shown states 
are closed channels. Note that there is no line-to-line relation between the states 
listed in each column. 
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Abstract. We analyse the formation of ultracold 7 Li 133 Cs molecules in the 
rovibrational ground state through photoassociation into the B 1 !! state, which 
has recently been reported [J. Deiglmayr et ah, Phys. Rev. Lett. 101, 133004 
(2008)]. Absolute rate constants for photoassociation at large detunings from 
the atomic asymptote are determined and are found to be surprisingly large. 
The photoassociation process is modeled using a full coupled-channel calculation 
for the continuum state, taking all relevant hyperfine states into account. The 
enhancement of the photoassociation rate is found to be caused by an "echo" of 
the triplet component in the singlet component of the scattering wave function at 
the inner turning point of the lowest triplet a 3 £ + potential. This perturbation 
can be ascribed to the existence of a broad Feshbach resonance at low scattering 
energies. Our results elucidate the important role of couplings in the scattering 
wave function for the formation of deeply bound ground state molecules via 
photoassociation. 
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1. Introduction 

The formation and manipulation of molecules at ultralow temperatures has undergone 
tremendous progress in the past years [TJ E]- The 'holy grail' of the field is the 
creation of a stable quantum-degenerate gas of molecules, in which all molecules 
populate the ground state of the system, comprising motional and internal degrees 
of freedom. Such a system provides exquisite possibilities for further control and may 
thus serve as the starting point to investigate future applications of ultracold molecules 
such as precision measurements of fundamental constants [3 and ultracold chemical 
reactions [I]. Dipolar molecules are of particular interest in this context, as the dipolar 
interaction can be exploited for the study of quantum many-body phenomena [S] [3] 
or the realization of various schemes for quantum information processing [TJ [8l [9] [10] . 
The most promising path towards the goal of a quantum-degenerate gas of molecules 
in the internal ground state consists in the association of ultracold atoms, either by 
magnetic or oscillatory electric fields. In the first case, called magnetoassociation, the 
molecules are formed through a magnetically induced Feshbach resonance coupling the 
free pair of atoms to a weakly bound molecular state. The transfer into the absolute 
internal ground state can then be accomplished by stimulated rapid adiabatic passage 
(STIRAP) in a combination of pulsed laser fields. This scheme has recently been 
successfully applied to form deeply bound molecules of Rb 2 [H], Cs 2 [HI H3], and 
KRb [14 . The other, closely related approach consists of the direct photoassociation 
of molecules out of an ultracold gas followed by spontaneous emission, which may be 
followed by a second bound-bound excitation step to transfer the molecules into the 
vibrational ground state. This method has lead to the formation of the vibrational 
ground state of K 2 [15], RbCs [16], Cs 2 [12], and LiCs [18]. Magnetoassociation 
combined with STIRAP has the great advantage of being a fully coherent process 
which preserves the phase-space density of the initial gas, while photoassociation 
followed by spontaneous emission leads to smaller phase-space densities but can be 
driven as a continuous process thus allowing for the steady accumulation of molecules, 
e.g. in an optical or static trap. Concerning ultracold dipolar gases we note that 
among the above listed ultracold molecules in the ground state the only dipolar ones 
are the heteronuclear molecules RbCs, KRb, and LiCs. They all posses significant 
dipole moments of 1.2 [19], 0.57 [14] , and 5.5 Debye [19] respectively. 

Here we focus on a recent experiment on the formation of ultracold bosonic 
7 Li 133 Cs ground state molecules with low vibrational and rotational quantum 
numbers [18[ 120) . The molecule formation consists in a particularly simple two-step 
photoassociation (PA) procedure comprising the laser excitation of the B 1 ]! state out 
of an ultracold mixture of lithium and cesium atoms. The vibrational state v'=4 of 
this electronically excited state has a particularly large Franck-Condon overlap with 
the vibrational level v"=0 of the ground state X 1 Y> + which leads to a measurable 
population of this lowest vibrational level by spontaneous emission. As only low 
angular momenta are involved due to the ultralow temperature, ultracold molecules 
in the rovibrational ground state could be detected. The rate limiting process of this 
PA scheme was the excitation of the tightly bound B 1 Ii(v' =A) level out of the dilute 
gas. Actually, the observation of PA of molecules in this state at short internuclear 
distance came as a surprise, as naive estimates delivered negligible formation rates. In 
this article, we investigate the PA process in more detail. We find that the PA rate is 
strongly enhanced by an increased amplitude of the scattering wave function at short 
internuclear distances. It is important to note that despite the presence of spin-orbit 



Influence of a Feshbach resonance on the photoassociation of LiCs 



3 



c 

LU 




4 



6 



8 



10 



12 



14 



R[A] 



Figure 1. Schematic view of the photoassociation process at short internuclear 
distances: the component of the scattering wave function related to the singlet 
potential X 1 E+ (solid black line) is coupled around the inner turning point of the 
triplet potential a 3 S+ to a level in the excited singlet state B 1 !! (black dashed 
line). A typical bound state wave function in the a' j S+ potential is also shown 
(gray line). 



coupling close to the atomic asymptote, for deeply bound levels the excited state has 
very pure singlet character. Therefore only the singlet component of the scattering 
wave function is relevant for the PA rate. Figure Q] depicts this situation: the Franck- 
Condon overlap between the scattering wave function in the singlet potential X 1 12 + 
and low lying levels in the B 1 !! potential vanishes, as the fast oscillations of the 
scattering wave function at short internuclear distances average to zero. However, 
experimentally we observe high PA rates into excited levels, which correspond to a 
transition around the inner turning point of the a 3 E + state. Therefore we conclude 
that the scattering wave function is locally perturbed by the presence of a bound 
level with strong triplet character, indicating the proximity of a Feshbach resonance. 
Indeed the importance of closed coupled channels for the scattering wave function 
close to a Feshbach resonance has recently been theoretically explored by some of us 
as a means to strongly enhance the PA formation rate of ultracold molecules |21) . 
The perturbation introduced by this Feshbach resonance leads to an "echo" of the 
triplet-like wave function on the singlet component of the scattering wave function, 
which increases the Franck-Condon overlap between the scattering wave function and 
tightly bound levels in the B 1 Ii potential. We present a model containing an accurate 
and complete description of the continuum scattering wave function and of the excited 
molecular wave function, which proves the validity of this picture and reproduces the 
experimental observations accurately. 

The paper is organized as follows: In section[2]we describe our experiment on the 
formation and detection of ultracold LiCs molecules in low vibrational and rotational 
states of the singlet ground state. Section [3] is then devoted to a quantitative analysis 
of measured line strengths and rate constants for PA into the B 1 !! state. In section^] 
we develop the theoretical model: A brief overview is given on how PA rates are 
computed based on wave functions of the excited levels. It is explained, how hyperfine 
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interactions in the scattering entrance channel couple singlet and triplet components 
and lead to an "echo" of triplet features in the singlet component of the scattering 
wave function. Numerical methods are introduced, and the results obtained by the 
coupled-channel calculation are discussed. Results for the coupled and uncoupled 
cases are compared with the experimental findings. Section [5] discusses the results of 
the paper. 

2. Formation and detection of ultracold LiCs molecules 

Details of the experimental setup for the formation and detection of ultracold LiCs 
molecules have already been described in detail in Refs. [T8j [20j [22] . Therefore, we will 
confine ourselves to the description of the main experimental features. 4 x 10 7 133 Cs 
atoms and 10 8 7 Li atoms are trapped in overlapped magneto-optical traps (MOT's) at 
densities of 3x 10 9 cm -3 and 10 10 cm -3 respectively. The MOT for cesium is realized as 
a dark spontaneous force optical trap [23 , in which the atoms are kept most of the time 
(typically 97%) in the dark lower hyperfine ground state leading to higher densities 
and to a reduction of inelastic collisions. In the conventional lithium MOT, most of 
the atoms populate the upper hyperfine ground state (>80%). Lithium and cesium 
atoms therefore collide mainly on the Li(2 2 Si/ 2 ,/=2)+Cs(6 2 Si/ 2 ,/=3) asymptote. 
We measure a cesium temperature of 250(50)/iK using time-of-flight expansion and 
deduce a lithium temperature of 600(150)/iK by fitting the line shape of a narrow 
PA resonance [20] to the model of reference [24]. In the center-of-mass frame, this 
corresponds to a mean collision energy of 580(80)/iK. 

In PA, a colliding pair of atoms absorbs a photon resonant to a transition into 
a bound excited molecular level [35] [55] . For this process we use up to 500 mW of 
light from a tunable Ti:Sa laser. The beam is collimated to a waist of 1.0 mm and 
continuously illuminates the two overlapped atom traps. The excited molecules decay 
within a few tens of nanoseconds either into bound ground state molecules or back into 
free pairs of atoms with additional kinetic energy [27]. For the detection of ground 
state molecules, we use a pulsed dye laser with a repetition rate of 20 Hz (typical 
pulse energy 8mJ, beam diameter ~5mm). Molecules in the ground state are ionized 
by resonant-enhanced multi-photon ionization (REMPI), in our case by two photons 
of the same color. The resulting LiCs + ions are then separated from a background 
of Cs + ions in a time-of-flight mass spectrometer and are finally detected using a 
micro-channel plate and single-ion-counting electronics. Figure[5](a) depicts the states 
involved in the formation and detection of ground state molecules. Exemplary PA 
resonances are shown in figure [2](b). 

We focus on the B 1 !! state of LiCs, which has been studied previously in a heat- 
pipe with high- resolution laser-induced fluorescence spectroscopy [28 . Due to spin- 
orbit interaction, this state is asymptotically correlated to Li(2 2 Si/2)+Cs(6 2 P3/2) after 
its coupling with neighboring triplet states. For convenience, we will refer to it as the 
B 1 Ii state. Additionally to the vibrational levels v' = — 25 identified in reference [25J, 
we found vibrational levels v' = 26 — 35, where the last level v' — 35 is bound only 
by a few GHz. For most lines, we observe rotational components J' — 1 and J' = 2, 
which are split into several hyperfine sub-lines. The typical overall width of these 
hyperfine structures ranges for J' = I resonances from 550 MHz for low vibrational 
levels to 900 MHz for the last bound levels, while the J' = 2 resonances are narrower, 
ranging from a width of 300 MHz for low v's to 600 MHz close to the asymptote. 

We detect PA resonances by ionozing ground state molecules produced after 
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Figure 2. (a) Photoassociation and detection scheme together with the 
relevant levels of LiCs: after PA into levels of the B 1 H state, molecules decay 
spontaneously into X S+ molecules. These molecules are ionized by two photons, 
where resonant intermediate states enhance the ionization rate, (b) Exemplary 
photoassociation resonances for deeply bound levels (u'=l) and levels close to the 
asymptote (u'=31) (detected at different REMPI wavelengths, see text). 



spontaneous decay from the photoassociated, electronically excited molecules. Three 
different frequencies were used for the ionization of these ground state molecules. 
Molecules produced by PA in low vibrational states v'<9 are ionized at 16999. 7 cm -1 , 
where three resonant transitions for the first step, X 1 £ + ,u"=0— KB 1 n,u / =14, 
X 1 S+,u"=l->- J B 1 n,u / =18, and ,v"=2-*B 1 TI,v'=23, are nearly degenerate |20] . 

Ground state molecules formed after PA via intermediate vibrational levels v'> 9 and 
v'< 15 are detected at 16859.4 cm" 1 , a REMPI resonance where a strong contribution 
from the transition X 1 E + ,v"=2->B 1 IL,v' = 19 was identified [20]. While the here 
identified ground state levels make up for an important part of the observed ion signal, 
molecules in other ground state levels might also be ionized by the given laser frequency 
via intermediate states other than B II. Molecules which have decayed from high lying 
vibrational levels v' > 15 are all ionized at 14692.7 cm -1 . The ionization process at this 
resonance has not yet been analysed. However, it is very likely that high lying levels 
in the electronic ground state X l Y, + are ionized at this energy through a broad band 
of levels in different electronic states between the asymptotes Li(2 2 S 1 /2)+Cs(6 2 P 3 / 2 ) 
and Li(2 2 S 1/2 )+Cs(5 2 D 3/2 ). 
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Figure 3. Photoassociation line strength (given as number of 

molecules X MHz/s/mW) and PA rates of all observed levels v',J'=l (left, 
dark colors) and J'=2 (right, light colors) in the B 1 !! state. The PA resonances 
are integrated and normalized to the PA laser power. Absolute rates are given 
for a reference PA intensity of 30W/cm 2 . Ionization wavelengths used for 
detection are 14692.7cm -1 (red), 16859.4cm -1 (green), and 16999.7cm -1 
(blue), respectively. 

3. Line strengths of the photoassociation resonances and formation rate 
constants 

3.1. Photoassociation resonance line strengths 

The strength of the PA resonances are shown in figure [3] for all observed v' , J'=l 
and J'=2 levels. The last bound level z/=35 was excluded from the analysis, as 
the cesium MOT is strongly perturbed by the PA laser at small detunings from the 
atomic asymptote. In order to reduce the influence of the varying hyperfine structure, 
the full PA spectrum is background-subtracted and integrated for each ro- vibrational 
level, therefore summing up the contributions from all hyperfine components. In 
order to directly compare the line strength of PA resonances measured under different 
experimental conditions, the integrated peak area is normalized to the intensity of 
both the PA laser and the ionization laser. At a strong line we confirmed for 
both the PA and the REMPI step that the molecule formation rate is indeed linear 
in both intensities. For the two-photon ionization, this indicates that the first 
resonant-enhanced transition is strongly saturated and therefore the probability for 
this transition does not depend on the intensity of the ionization laser above a certain, 
not determined threshold. The PA laser intensities were chosen for each resonance in 
such a way that the maximum count rate does not exceed 0.5 ions/pulse to avoid 
saturation of the single ion detection. 

One uncertainty remains, which limits the comparability between different PA 
resonances. As discussed in Sect. [2j we do not detect the excited molecules directly 
but only the ground state molecules formed after spontaneous decay. As the Frank- 
Condon factors between excited state and ground state levels vary with the PA level 
v', also the distribution of populated ground state levels after spontaneous decay 
depends on the PA level [20]. These differences in the populated levels can lead to 
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variations in the overall ionization efficiency for molecules produced at different PA 
resonances. It is difficult to estimate the influence of this varying detection efficiency 
on the measured PA strengths, as many off-resonant and nonlinear processes play a 
role in the REMPI ionization scheme. We observe however that the calculated ground 
state distributions vary only slowly with the excited level v'. Also from measurements 
of the same PA line at different ionization resonances, we deduce that this variation 
does not exceed a factor of ten over neighboring vibrational levels. Nonetheless this 
constitutes one of the main uncertainties in the reported PA rate constants. For levels 
close to the asymptote we also compare the molecule formation rate deduced from the 
ion signal with measurements of trap loss. We find fair agreement within one order 
of magnitude, which is acceptable in view of the large systematic uncertainty on the 
efficiency of the ion detection. 

3.2. Absolute molecule formation rates 

In order to relate the PA line strengths to the formation rate of molecules, we further 
analyse the detection efficiency for ground state molecules. As the formed molecules 
are not trapped, they leave the ionization region some time after their production due 
to their thermal velocity of roughly 300 fiK (calculated from the velocities of the atomic 
constituents) and acceleration due to gravity. Taking the size and the alignment of 
the ionization beam (roughly one beam diameter below the trapped atom clouds) 
into account, we estimate a geometric overlap factor of 40%. Therefore, only 40% of 
the formed molecules have a chance to be ionized. Assuming that the first bound- 
bound transition during the ionization of ground state molecules is on resonance and 
thus fully saturated, we use a typical ionization cross section of 10 _18 cm 2 [29] to 
approximate the ionization probability p lon for pulse energies Erempi of up to few 
tens of mJ by p lon ~ 8x I _3 Erempi where Erempi is measured in mJ. Finally the 
detector efficiency is around 20% [30], leading to an overall efficiency of the detection 
setup on the order of 10~ 3 for typical experimental parameters. By averaging over 
many cycles (typically around 100), even very weak PA resonances can be detected. 

The rate constant ifpA for the PA process is calculated by first linearly scaling 
the number of produced molecules per second to a reference PA intensity (chosen to be 
30W/cm 2 ) yielding a molecular formation rate k. In our experimental geometry the 
diameter of the PA laser beam and of the lithium MOT are both larger than the size of 
the cesium MOT. The rate coefficient Kpa is then simply given by ifpA = k/n^i/Ncs 
with the peak lithium density n^i and the cesium particle number Nc s - The measured 
rate coefficients are shown in figure [3] (right axis). For the last bound states we observe 
PA rate constants on the order of 4x 10~ n cm 3 /s. These are of the same order of 
magnitude as the rate constants observed in the PA of Cs2 [Ml [32] and KRb [33] , as 
was predicted by Azizi et al. |34j . 

The PA rate decreases towards lower vibrational states, which are excited at 
shorter internuclear distance corresponding to lower pair density. Below v' =25 we 
observe a change towards a distinct oscillatory pattern in the rate constants. While 
some vibrational levels show rate constants varying by no more than one order of 
magnitude, adjacent levels are below the detection threshold of our setup. For very 
deeply bound levels v'<6 the PA rate appears to increase again, e.g. for the v'=l,J'=l 
level at a detuning of more than 1.400 cm -1 we measure a PA rate coefficient on the 
order of 4xl0~ 13 cm 3 /s. This is unexpectedly high, as a single channel scattering 
calculation predicts a rate constant below 10~ 19 cm 3 /s due to the low Franck-Condon 
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overlap between the rapidly oscillating singlet scattering wave function and the excited 
state wave function at short internuclear distances. In the following, we describe the 
theoretical framework employed to explain the unexpectedly large formation rate in 
the singlet electronic state. 



4. Model of photoassociation into the B 1 !! electronic state 

4--1- Theoretical photoassociation rates 

A PA rate Kp A = {v rc \crp A } from an initial scattering continuum state to an excited 
vibrational level v' is obtained using the cross section cr v p A of two atoms colliding 
in the presence of a laser held, and averaged over the relative velocity distribution 
Urei- For a thermal atomic gas, a Maxwell-Boltzmann distribution characterized by a 
temperature T is appropriate, and assuming a PA laser beam of negligible linewidth 
and resonant for the transition to v' , we have |35j : 

K$AT,I) = -r-pr — / de e- e ' kBT , (1) 

where Q T = {2irfj,k B T/h 2 ) 3/2 , h and k B are the Planck and Boltzmann constant, 
respectively, e = piv 2 e j2 is the relative kinetic energy of the colliding pair of atoms of 
reduced mass /it, and j v i is the natural linewidth of the photoassociated level. 

Using Fermi's Golden Rule, we obtain an expression for the stimulated emission 
width 7 S in terms of the initial scattering state ^ Cj ; at energy e, and the excited 
vibrational state 4> V ',j' 

ls = —\(<f> v ,,r\D(R)\* e ,e}\ 2 . (2) 

Here, I labels the partial wave, I is the intensity of the PA laser, and eo and c arc 
the vacuum permittivity and speed of light, respectively. We will restrict ourselves 
to s-wave scattering in the entrance channel, and thus set £—0. This is well justified, 
as the measured relative temperature in the experiment of 580 fiK is well below the 
p-wave barrier height of 1.6 mK. 

We obtain the scattering wave function in the entrance channel ^ c ./ by solving 
the Hamiltonian of two colliding atoms with hyperfine interactions, which mix the 
singlet A X S + and triplet a 3 E + states. In addition to this scattering wave function 
and the wavefunction of the excited level <fi v >,j> (see descriptions in the next sections), 
we need the dipole transition moment D(R) to calculate j s . Its i?-dependent form 
is calculated using the theoretical method and parameters described in reference [19] . 
The value of the dipole moment does not vary by more than 10% around a mean value 
of 4.0 a.u. over the relevant range of distances. 



4-2. Excited state vibrational wave functions 

PA is performed to the excited B 1 !! molecular potential. High resolution spectroscopy 
of this potential has been published in reference [28] . We use a potential fitted to the 
levels observed in this work combined with the additional levels observed by us |36j . 
However, we note that the spin-orbit coupling of the B 1 !!, b 3 H, and c 3 E + states 
becomes important near the dissociation limit and cannot be neglected. The Hund's 
case (a) is therefore no longer valid and the B 1 !! state should be treated as f2=l in 
Hund's case (c). To estimate the strength of the singlet-triplet mixing we used a simple 
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Figure 4. Projection of the Q=l (B 1 !!) wave functions onto the singlet 
component B 1 !! (light blue) and the triplet components fe 3 II and c 3 S+ (dark 
red). At low v's the state can therefore be identified as pure B 1 ]!, while for 
i/~20 and higher the triplet character, induced by spin-orbit coupling, becomes 
important. 



3- channel model [37], coupling ab-initio curves for B 1 !!, 6 3 LT, and c 3 £ + (calculated 
as described in reference |19j ) with constants proportional to the atomic finestructure 
interaction of 554.04 cm -1 [3H1 f° r the first asymptote of cesium. Wave functions 
and projections on the three channels where calculated using the mapped Fourier grid 
method 0D]. 

As can be seen in figure 0] the components of the £1=1 (B 1 !!) wave function 
related to the coupled triplet states become sizable above v' — 20. For the highest 
lying levels, the summed amplitude of these components reaches nearly 50%. This 
value can be taken as a first estimate for the influence of the excited state mixing on 
the PA line strength: the triplet component of the excited state can now be populated 
from the triplet component of the scattering state. In general (i.e. as long as there are 
no strong local perturbations of the wave functions), the excited triplet components 
will only decay into triplet ground state molecules or free pairs of atoms, so they will 
not contribute to the production of deeply bound singlet molecules. In such a case 
one would expect for high lying levels to observe PA rates, which are reduced by the 
amount of the triplet projection shown in figureS]when compared to PA rates derived 
from the here presented theory. However, since the focus of this work is on the lowest 
vibrational levels (below v' = 20), we treated the excited potential as a pure B 1 !! 
state in our calculation. 

4- 3. Scattering wave function of the ground state 

If we label sj and ij the electronic and nuclear spin for the atom j, respectively, 
then the total spin fj of atom j (with projection mj^ is given by fj — sj + ij. The 
total angular momentum without rotation is F = f± + /a. The two-body hypcrfinc 
Hamiltonian for atoms of relative momentum p (and reduced mass fi) can be written 
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as 

„2 2 n U) 

where y c is the potential energy operator, and a$ is the hyperfine constant for atom 
j. No spin- rotation coupling is introduced here, which means that the rotational 
angular momentum I is fixed. V c can be expressed in terms of the singlet Vo(R) and 
triplet Vi(R) molecular potentials using the projection operators P° and P 1 : 

V c = V {R)P° + VyifyP 1 . (4) 

As it is the case for the excited B 1 Ii state, the ground electronic singlet and triplet 
potentials are known from high resolution spectroscopy [H]. We used the potential 
curves as given in the reference, with a cubic spline interpolation of the given pointwisc 
representation connected to the given long-range expansion. 

The hyperfine Hamiltonian in Eq.Q only couples hyperfine channels with the 
same total angular momentum projection rap = m f\ + 171 f 2 - F° r 7 Li with i\ = 3/2 
and 133 Cs with i 2 = 7/2 (both with si = s 2 = 1/2), mp takes values from -6 to +6, 
leading to a total degeneracy of the atom pair of 128. However, atoms mainly collide 
on the Li(2 2 Si/ 2 ,/ = 2)+Cs(6 2 Si/ 2 , / = 3) asymptote, and thus only 35 degenerate 
entrance channels have to be considered. 

The hyperfine Hamiltonian ([3]) is expressed and diagonalized in the atomic basis 

\h,m 1 ;f 2 ,m 2 ) = |/i,mi)Li ® |/2,m 2 )c s , (5) 
and the total collisional wave function is then expressed as 

N 

|*M> = E MR){\fum x ) ® |/a, m 2 )} a , (6) 

a=l 

where a labels a particular channel, and N is the number of coupled hyperfine channels 
which depends on the value of mp. Here, -0 Q (i?) is the radial wave function of channel 
a. 

For the calculation of the Franck- Condon factors, only the singlet component of 
the total collisional wave function is relevant. A rotation is thus performed to express 
the wave function in the molecular basis \S, ms, I, mi) with S = s*i+s*2 and I = i\ +i 2 - 
We then have 

N 

|*M> = J2MR){\S,m s ,I,m I )}p , (7) 

13=1 

where (f)p(R) is the wave function in this basis (with TV channels (3). We can go from 
one basis to another basis using angular momentum algebra. 

4-4 ■ Results for wif = 5 scattering state 

For simplicity we concentrate here on the particular case of mp = 5, corresponding 
to the coupling of the four channels listed in Table [1] by the hyperfine interaction. 
These four states represent approximately l/10th of a total of 35 degenerate channels 
of the entrance channel. In order to compute the dipole transition matrix elements 
(</v,j'|£ ) (-ft)l , I't,i) appearing in Eq.@ for 7 S , the wave functions are needed: the 
Schrodinger equation is solved with the Mapped Fourier Grid Hamiltonian (MFGH) 
method 40, 21 . All calculations are performed at a collision energy of 500 jjK: No 




Figure 5. Components of the energy- normalized radial wave function for a Li-Cs 
pair colliding at 500 fiK, in the atomic basis (left columns) and in the molecular 
basis (right column). Under the given experimental conditions, the upper- most 
atomic state /i=2,m/i=2,/2=3,m/2 = 3 is the only open entrance channel, the 
three other shown states are closed channels. Note that there is no line-to-line 
relation between the states listed in each column. 



large variations of the results are expected for slightly varying collision energies, so the 
experimental collision energy is very well approximated by this value. The components 
of the initial radial wave function are drawn in figure [5] in the atomic (left column) and 
molecular (right column) basis. It is striking to see in both columns that, in contrast 





Atomic basis 


Molecular basis 


rriF = 5 


fx 


= 2, to/1 = 2, f 2 = 3, 


m/2 


= 3 


S = 


0, ms = 0,7 = 5, mi 


= 5 




fl 


= 2, to/1 = 2, h = 4, 


TO/2 


= 3 


s = 


1, ms = 0,7 = 5, to/ 


= 5 




fl 


= 1,TO/1 = l,/ 2 = 4, 


TO/2 


= 4 


s = 


1,TOS = 1,1 = 4, TO/ 


= 4 




fl 


= 2, to/1 = l,/ 2 = 4, 


TO/2 


= 4 


s = 


1, TOs = 1,7 = 5, TO/ 


= 4 



Table 1. Values for the quantum numbers of the collisional channels of the 
mp = 5 subspace for the atomic basis |/i , m/i , /2, m f2) an d the molecular basis 
\S, ms, I, mi). Note that there is no line-to-line relation between the states listed 
in each column. 
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Figure 6. Black lines: Singlet projection of the total wave function in the 
coupled-channel model (upper panel), and in a single channel calculation with 
the singlet potential (lower pannel). Probability densities for the v'=l level of 
the B l H state are also displayed (red dashed lines), illustrating the increased 
overlap due to the enhancement of the ground state wave function at the inner 
turning point of the triplet potential. 



with the generally accepted picture of such a coupling case, the components exhibit 
clear irregular features resulting from the hyperfhic coupling even at short internuclear 
distances, i.e. well inside the region where the atomic states decouple into triplet and 
singlet states. By looking at the molecular basis decomposition (right column), it is 
clear that the singlet component (upper panel) is so strongly coupled to the triplet 
component {e.g. third panel from the top) that they contaminate each other with their 
amplitude variation. In particular, the amplitude of the singlet component is strongly 
enhanced in the range of the triplet potential, i.e. from R = 7 a.u. towards large 
distances, compared to the uncoupled singlet wave function (see figureEJ. This "echo" 
of the triplet component appearing on the singlet projection of the wave function 
indicates that bound levels of the closed channels are energetically close to the open 
entrance channel. Only this situation explains the strong local perturbation of the 
scattering wave function. Such a hyperfine-mediated coupling between a bound level 
and a zero-energy continuum state is generally referred to as a Feshbach resonance. 

In order to analyse how close the modeled situation is to the Feshbach resonance, 
we look at the squared amplitude of the collisional wave function on the closed channels 
as a function of the scattering energy. With the total state \^ e ,e) of the coupled system 
at energy e being expressed as \^ e ,e) — 2 a =i V^,a(-ff) \oi), where a labels the atomic 
hyperfme states, we define the amplitude P a (E) of a given hyperfine state at energy e 
as P a (E) — J Q \ip eyCC (R)\ 2 dR. This definition is only meaningful for closed channels, 
where P a =0 corresponds to an unpopulated channel and P Q =1 indicates, that \^ e ,l) 
is a pure bound state in channel a. We plot the results for the three closed channels 
in figure We see that indeed a resonance is present in the upper panel around 
6 x 10~ 7 a.u. (for t = 0) and around 8 x 10~ 7 a.u. (for I = 2), as maxima in the 
amplitude of the closed-channel wave functions. Traces of these resonances are still 
present close to zero collision energy, where the experiment takes place. Therefore, the 
enhancement of the amplitude of the total scattering wave function is induced by the 
increased mixing arising from the presence of a broad Feshbach resonance at higher 
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Figure 7. Amplitude of the mp=5 initial wave function on the three closed 
channels as a function of collisional energy, for £=0 (upper panel) and 1=2 (middle 
panel). The inset in the upper graph shows an enlargement of the energy range 
for £=0 including the collision energy of 500 /iK (1.6xl0~ 9 a.u.), used in all other 
calculations. The amplitude is given as P a (E), the weight of the component a 
associated to one hyperfine channel in the total wave function. Details are given 
in the text. Lower panel: elastic cross section for the same cases mp=5, £=0,2. 
The positions of Feshbach resonances are indicated by vertical dashed lines. 



energy. In particular the component of the f\ — \,mf 1 = l,/2 = 4,m/ 2 =4 state 
dominates at low energies. Note that the exact position of the resonance might not be 
too reliable due to uncertainties especially in the used triplet ground state potential, 
as discussed by the authors of reference |41) . However as we will show later, even the 
very weak mixing introduced by the Feshbach resonance in this calculation leads to a 
convincing reproduction of the experimentally observed rates, so the precise position 
of the resonance seems of minor relevance. 

We also calculate elastic cross sections for all mp subsystems using the log- 
derivative algorithm of Johnson [32]. This is another way to locate and identify 
possible resonances either due to the shape of potentials (shape resonances), or due 
to the coupling with close-by bound states (Feshbach resonances) . In figure [7] (lowest 
panel) we show the energy-dependent elastic cross section for mF—5. Also here, 
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superimposed on the regular nodal structure of the elastic cross section as one would 
expect it for an uncoupled case, the trace of the above mentioned resonances is visible. 
It is important to note that this Feshbach resonance is actually present at about the 
same energy in the elastic cross-sections for all values oi mp. This suggests that the 
singlet component of the wave function will be perturbed for all mp in a similar way 
as in the presented tuf = 5 case. Therefore the enhancement of the amplitude in 
the singlet channel should be observable for all values of mj? in a similar way and we 
expect that the PA rate averaged over all experimentally possible values for tuf will 
not strongly deviate from the value derived here for mp=5- 

4-5. Comparison with experimental photoassociation rate constants 

The PA rate is proportional to the transition dipole moment matrix element between 
the singlet component of the initial continuum state with the vibrational levels of the 
B state. The calculated rates assuming identical contributions from all mp subsystems 
is represented in figure [3] for both coupled and uncoupled initial wave functions. As 
expected, both cases yield similar rates for high vibrational states, induced by a similar 
long-range part of the wave function. For low levels of the B state, however, the 
PA rate is strongly enhanced in the coupled-channel picture due to the "echo" of 
the triplet-like wave function on the singlet component of the wave function. It is 
evident that the oscillatory pattern of the experimental PA rate is well reproduced 
by the coupled-channel model. This oscillatory pattern arises from the variation of 
the overlap between the initial scattering wave function and the wave function of the 
excited level as the outer lobe of the latter "scans" with changing v' over the nodal 
structure of the former. Fair agreement between the absolute rates from experiment 
and theory is also found. Even though the measured values suffer from a systematic 
uncertainty of roughly one order of magnitude, it is obvious that the experimental 
observations are not compatible with the conventional picture of uncoupled singlet 
and triplet states at short internuclear distances. 

5. Conclusion 

By analyzing photoassociation rates of LiCs molecules in the B 1 !! state we find that 
the formation rate of ground state molecules is governed by a perturbation of the 
ground state scattering wave function. The perturbation is caused by the existence 
of a Feshbach resonance at low energies, which enhances the PA rate drastically. 
Specifically it was shown that the singlet component of the scattering wave function, 
which is responsible for the efficient transfer into the X 1 Y 1 + ground state state through 
photoassociation and subsequent spontaneous emission, contains a strong "echo" of 
the triplet wave function, which allows for an efficient transfer of population from 
free pairs into bound molecules. As the photoassociation rate into the strongly bound 
levels of the excited B 1 Ii state is the rate limiting process for the formation of deeply 
bound ground state molecules, the perturbation leads to an increase of the formation 
rate by several orders of magnitude when compared to the uncoupled case. 

To which extent our observations are specific to the LiCs case remains to be 
investigated in future experiments, but we speculate that a similar enhancement can 
be found in all other alkali systems. This is supported by the fact that the observed 
strong enhancement of the molecular formation rate is induced by the presence of 
a Feshbach resonance roughly 200mK (or 4 GHz) above the experimental collision 
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Figure 8. Experimental photoassociation rate constants (same data as figure(3]l 
for the J'=l lines (filled squares), compared to the predicted values obtained for 
mp = 5, in the coupled-channel picture (solid line) and in the uncoupled case 
(dashed line). Vibrational levels, which could not be detected experimentally, 
are marked by arrows indicating the detection threshold of the experiment. In 
the model calculation, a PA laser intensity of 30W/cm 2 and a collision energy of 
500 /jK was used, corresponding to the experimental parameters. Note that the 
systematic error on the experimental values is about one order of magnitude, as 
described in Sect. 13.21 



energy, a situation which is not uncommon in alkali mixtures. We note that this 
kind of enhancement, or R-transfer, through coupling in the scattering wave function 
is complementary to other proposed and realized schemes for the transfer of weakly 
bound triplet molecules into ground state singlet molecules [43]. In these schemes, 
the "transfer" from triplet to singlet character occurs via spin-orbit mixing in the 
electronically excited state. In RbCs and KRb, levels in the spin-orbit coupled states 
c 3 E + and B 1 !! were used to transfer weakly bound triplet molecules into absolute 
ground state molecules with pulsed laser fields [THl [H] , while in Cs2 the coupled state 
0„ (b 3 II„-A 1 S+) was successfully employed in a STIRAP step towards deeply bound 
molecules [12]. We suggest that also the wave function of most so called Feshbach 
molecules might exhibit an enhanced amplitude in the singlet channel at the inner 
turning point of the triplet potential. This would facilitate the transfer of these 
molecules into absolute ground state molecules via suitable excited singlet states. 

The observed strong enhancement of the PA rates into the rotational component 
J'=2 as compared to the component J'=l for low v's is still an open question which 
requires further investigations. A possible explanation would be an enhancement due 
to a d-wave shape resonance [18j . while no indication for such a resonance was found 
in the model presented here. Another reason for the observed ratio between rotational 
states could be rotational coupling with higher partial waves, which was neglected in 
our model. In fact, the strongly differing widths for J'—l and J'=2 lines could indeed 
be an indication of the necessity to take into account such rotational coupling at short 
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distances. 
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